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ABSTRACT 

We study of the effect of turbulence on heat transfer within magnetized plasmas for energy injection 
velocities both larger and smaller that the Alfven speed. We find that in the latter regime the heat 
transfer is partially suppressed, while in the former regime the effects of turbulence depend on the 
intensity of driving. In fact, the scale I a at which the turbulent velocity is equal the Alfven velocity 
is a new important parameter. When the electron mean free path A is larger than I a, the stronger 
the the turbulence, the lower thermal conductivity by electrons is. The turbulent motions, however, 
induces their own advective heat transport, which, for the parameters of intracluster medium (ICM) 
provides effective heat diffusivity that exceeds the classical Spitzer value. 
Subject headings: turbulence - ISM: general - galaxies: clusters: general - MHD 



1. ASTROPHYSICAL MOTIVATION 

Heat transfer in turbulent magnetized plasma is an im- 
portant astrophysical problem which is relevant to the 
wide variety of circumstancies from mixing layers in the 
Local Bubble (see Smith & Cox 2001) and Milky way 
(Begelman & Fabian 1990) to cooling flows in intraclus- 
ter medium (ICM) (Fabian 1994). The latter problem 
has been subjected to particular scrutiny as observations 
do not support the evidence for the cool gas (see Fabian 
et al. 2001). This is suggestive of the existence of heating 
that replenishes the energy lost via X-ray emission. Heat 
transfer from hot outer regions is an important process 
to consider in this context. 

It is well known that magnetic fields can suppress ther- 
mal conduction perpendicular to their direction. The is- 
sue of heat transfer in realistic turbulent magnetic fields 
has been long debated (see Bakunin 2005 and references 
therein). An influential paper by Narayan & Medvedev 
(2001, henceforth NM01) obtained estimates of thermal 
conductivity by electrons using the Goldreich-Shidhar 
(1995, henceforth GS95) model of MHD turbulence with 
the velocity Vl at the energy injection scale L that is 
equal to the Alfven velocity Va, i.e. the turbulence with 
the Alfven Mach number M A = (V L /V A ) = 1. This is 
rather restrictive, as in the ICM Ma > 1 (see §4), while 
in other astrophysical situations Ma < 1- Below we dis- 
cuss turbulence for both Ma > 1 and Ma < 1 and com- 
pare the heat transfer by electrons to that by turbulent 
fluid motions. 

2. THERMAL CONDUCTIVITY: STATIC MAGNETIC FIELD 

2.1. Basics of heat transfer in magnetized plasma 

Following NM01, we initially disregard the dynamics of 
fluid motions on heat transfer, i.e. consider thermal con- 
ductivity induced by electrons moving along static mag- 
netic fields. Magnetized turbulence in the GS95 model 
is anisotropic with eddies elongated along (henceforth 
denoted by ||) the direction of local magnetic field. Con- 
sider isotropic injection of energy at the outer scale L 
and dissipation at the scale lj_ t mim where _L denotes the 
direction of perpendicular to the local magnetic field. 
NM01 observed that the separations of magnetic field 
lines for tq < l± m i n are mostly influenced by the mo- 



tions at the scale l± im i n , which results in Lyapunov-type 
growth: ~ tq exp(Z/iii m in)- This growth is similar to 
that obtained in earlier models with a single scale of tur- 
bulent motions (Rechester & Rosenbluth 1978, Chandran 
& Cowley 1998). This is not surprising as the largest 
shear that causes field line divergence is provided by the 
marginally damped motions at the scale around l± im i n - 
In NM01 ro is associated with the size of the cloud of 
electrons of the electron Larmor radius rLar,eiectr- They 
find that the electrons should travel over the distance 



jrr 



M l ±,min/rLar,electr) 



(1) 



to get separated by l± t min- 
Within the single-scale model which formally corre- 



sponds to L = I 



| ,min 



l±,min the scale Lrh is called 



Rechester-Rosenbluth distance. For the ICM parameters 
the logarithmic factor in Eq. Q is of the order of 30, and 
this causes 30 times decrease of thermal conductivity for 
the single-scale models 1 . In the multi-scale models with 
a limited (e.g. a few decades) inertial range the logarith- 
mic factor stays of the same order but it does not affect 
the thermal conductivity, provided that L ^> ln m i n . In- 
deed, for the electrons to diffuse isotropically they should 
spread from rLar,electr to L. The GS95 model of turbu- 
lence operates with field lines that are sufficiently stiff, 
i.e. the deviation of the field lines from their original 
direction is of the order unity at scale L and less for 
smaller scales. Therefore to get separated from the ini- 
tial distance of l±, m in to a distance L (see Eq. with 
Ma — 1), at which the motions get uncorrelated the 
electron should diffuse the distance slightly larger (as 
field lines are not straight) than \[2L (NM01, also see 
§2.3.), which is much larger than the extra travel dis- 
tance ~ 30Zii m j n . Explicit calculations in NM01 support 
this intuitive picture. 

2.2. Heat Transfer for Ma > 1 
Turbulence with Ma > 1 evolves along hydrodynamic 
isotropic Kolmogorov cascade, i.e. Vi ~ Vl(1/ X) 1 / 3 over 

1 For the single-scale model Lrr ~ 30L and the diffusion over 
distance A takes Lfm/L steps, i.e. A 2 ~ LfmL, which decreases 
the corresponding diffusivity coefficient n e i ec tr, single ~ A 2 /<5t by 
the factor of 30. 
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the range of scales [L, I a], where 

Ia « WaIVl? = LM A \ (2) 

is the scale at which the magnetic field gets dynami- 
cally important, i.e. Vi — V A . This scale plays the 
role of the injection scale for the GS95 turbulence, i.e. 
Vi ~ VkC/x/^i) 1 ) with eddies at scales less than l A get- 
ing elongated in the direction of the local magnetic field. 
The corresponding anisotropy can be characterized by 
the relation between the semi-major axes of the eddies 

l\\ ~ L{l±/L) 2 'Hl A \ M A > 1, (3) 

where |j and _L are related to the direction of the local 
magnetic field. In other words, for Ma > 1, the tur- 
bulence is still isotropic at the scales larger to I a, but 
develops (I^/Ia) 1 ^ anisotropy for I < Ia- 

For electron mean free path A > i^, electrons stream 
freely over the distance of I a- For electrons at dis- 
tance l±, m in to get separated by L the required travel 
is the random walk with the step I a, i.e. the mean- 
squared displacement of a thermal electron till it enters 
an independent large-scale eddy A 2 ~ 1\(LJIa), where 
L/Ia is the number of steps. These steps require time 
St ~ (L/l A )lA/Civ e i ectr , where v e iectr is electron ther- 
mal velocity and the coefficient C\ — 1/3 accounts for 
ID character of motion along magnetic field lines. Thus 
the electron diffusivity coefficient is 

K e l ec t r = A 2 / 'St f» ( l/3)l A V e lectr , I A < A, (4) 

which for Ia <C A constitutes a substantial reduction 
of conductivity compared to its Spitzcr (unmagnetized) 
value Kspitzer = ^v e i ectr - We assumed in Eq. (@J that 
L » 3(% min (see §2.1). 

For A < I a < L, K e i ec tr ~ ^fi^spitzer as both the L RR 
and the additional distance for electron to diffuse because 
of magnetic field being stiff at scales less than Ia are 
negligible compared to L. For Ia — > L, when magnetic 
field has rigidity up to the scale L, it gets around 1/5 of 
the Spitzer value according to NM01. 

Note, that even dynamically unimportant magnetic 
fields do influence heat conductivity over short time 
intervals. For instance, over time interval less than 
l\jC\K spitzer the diffusion happens along stiff magnetic 
field lines and the difference between parallel and perpen- 
dicular diffusivities is large 2 . This allows the transient 
existence of sharp small-scale temperature gradients. 

2.3. Heat Transfer for Ma < 1 

It is intuitively clear that for Ma < 1 turbulence 
should be anisotropic from the injection scale L. In 
fact, at large scales the turbulence is expected to be 
weal? (see Lazarian & Vishniac 1999, henceforth LV99). 
Weak turbulence is characterized by wavepackets that 

2 The relation between the mean squared displacements per- 
pendicular to magnetic field (y 2 ) and the displacements x along 
magnetic field for x < Ia can be obtained through the diffusion 
equation approach in §2.3 and Eq. JJJ. This gives (j/ 2 ) 1 / 2 ~ 



3 The terms "weak" and "strong" turbulence are accepted in 
the literature, but can be confusing. As we discuss later at smaller 
scales at which the turbulent velocities decrease the turbulence 
becomes strong. The formal theory of weak turbulence is given in 
Galtier et al. (2000). 



do not change their iy, but develop structures perpen- 
dicular to magnetic field, i.e. decrease l± . This cannot 
proceed indefinitely, however. At some small scale the 
GS95 condition of critical balance, i.e. 1\\/Va ~ l±/Vi, 
becomes satisfied. This perpendicular scale Itrans can 
be obtained substituting the scaling of weak turbulence 
(see LV99) Vi ~ VlH^/L) 1 / 2 into the critical balance 
condition. This provides Itrans ~ LM\ and the corre- 
sponding velocity Vtrans ~ VlMa- For scales less than 
Itrans the turbulence is strong and it follows the scalings 
of the GS95-type, i.e. Vi ~ V L {L /l^ 1 ^ M X J 3 and 

1,1 ~ L(l ± /L) 2 /Hl A 4/3 , M A < 1. (5) 

For Ma < 1, magnetic field wandering in the direction 
perpendicular to the mean magnetic field (along y-axis) 
can be described by d{y 2 )/dx ~ (y 2 )/l\\ (LV99), where 4 
in is expressed by Eq. J5J and one can associate l± with 
2{y 2 ) 

tf) 1 ' 2 ~ ^m M A' l±< hrans (6) 

For weak turbulence d(y 2 )/dx ~ LM A (LV99) and thus 

( W a )Va„il/V/ 2 M2, l±> Itrans. (7) 

Eq. © differs by the factor M\ from that in NM01, 
which reflects the gradual suppression of thermal conduc- 
tivity perpendicular to the mean magnetic field as the 
magnetic field gets stronger. Physically this means that 
for Ma < 1 the magnetic field fluctuates around the well- 
defined mean direction. Therefore the thermal conduc- 
tion gets anisotropic with the coefficient of thermal con- 
duction parallel to the mean field n\\ } eiectr ~ ^/^ K spitzer 
being larger than K±, e iectr for the thermal conductivity 
in the perpendicular direction. 

Consider the coefficient k± e iectr 

for M A < 1. As 
NM01 showed, electrons become uncorrelated if they are 
displaced over the distance L in the direction perpendic- 
ular to magnetic field. To do this, an electron has first 
to travel Lrr (see Eq. Q), where Eq. JSJ relates In m i n 
and lx,min- Similar to the case in §2.1, for L 3> 30l\\ j7nin , 
the additional travel arising from the logarithmic fac- 
tor is negligible compared to the overall diffusion dis- 
tance L. At larger scales electron has to diffuse ~ L 
in the direction parallel to magnetic field to cover the 
distance of LM\ in the direction perpendicular to mag- 
netic field direction. Therefore the separation of elec- 
trons over the turbulence driving scale L perpendicular 
to the magnetic field direction requires L/LM\ = M A 2 
random steps. If A <^ L the diffusion over L requires 
time St ~ L 2 /(M A D\\), where D\\ is the diffusion coeffi- 
cient which is w e / ec t r A/3. As a result 

K±.electr = L 2 / St M l/3v e l ectr \M%, M A < 1, (8) 

where we disregarded the distance to travel in the direc- 
tion perpendicular mean magnetic field, i.e. L, compared 
to the distance to travel parallel to magnetic field, i.e. 
LM A 2 . For M A of the order of unity this is not accurate 
and one should account for the actual 3D displacement 
(see NM01 and §2.1). 

4 The fact that one gets ln m i n in Eq. Q is related to the 
presence of this scale in this diffusion equation. 
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Fig. 1. — Sonic Mach number M s is ploted against the Alfven 
Mach number M^. The heat transport is dominated by the dy- 
namics of turbulent eddies is above the curve and by thermal con- 
ductivity of electrons is below the curve. Here A is the mean free 
path of the electron, L is the driving scale, and a = (rris/mp) 1 ^ 2 , 
/3 4. The panel in the right upper coner of the figure illustrates 
heat transport for the parameters for a cool core Hydra cluster 
using data from EV06 (point "CC"), "I" corresponds to the illus- 
trative model in EVP05. 



3. FLUID VERSUS ELECTRON MOTIONS 

Turbulent motions themselves can advectively trans- 
port heat. In Cho et al. (2003) we dealt with the turbu- 
lence with Ma ~ 1 and estimated 

Kdynamic ~ Cdyn^-^Li Ma > 1, (9) 

where Cdyn ~ 0(1) is a constant, which for hydro tur- 
bulence is around 1/3 (Lesieur 1990). For fully ionized 
non-degenerate plasma we assume Cdyn ~ 2/3 to ac- 
count for the advective heat transport by both protons 
and electrons 5 . Thus eq. © covers the cases of both 
Ma > 1 up to Ma ~ 1. For Ma < 1 one can esti- 
mate Kdynamic ~ cPuj, where d is the random walk of 
the field line over the wave period ~ u> . As the weak 
turbulence at scale L evolves over time r ~ w -1 , 
(y 2 ) is the result of the random walk with a step d, i.e. 
(y 2 ) ~ (ruj)d 2 . According to eq.® and 0, the field line 
is displaced over time r by (y 2 ) ~ LM a Vat. Combin- 
ing the two one gets d 2 ~ LM^Vlw -1 , which provides 
K dynamic ~ C dyn^-V^M^ which is similar to the diffu- 
sivity arising from strong turbulence at scales less than 

hransi i-C. ^dynamic ~ CdynhransVtrans ■ The total diffu- 

sivity is the sum of the two, i.e. for plasma 

Kdynamic » (0/3)LV L M\, M A < 1, (10) 
where j3 w 4. 

The schematic of the parameter space for n e icctr < 
^dynamic is shown in Fig 1, where the the Mach num- 
ber M s and the Alfven Mach number Ma are the 
variables. For Ma < 1, the ratio of thermal con- 
ductivities arising from fluid and electron motions is 
^dyn amic/ Keiectr ~ (3aM s M A {L/X) (see Eqs. © and 
I|l6|0- the square root of the ratio of the electron to pro- 
ton mass a = (nie/mp) 1 ^ 2 , which provides the separa- 
tion line between the two regions in Fig. 1, (3aM a ~ 

5 This gets clear if one uses the heat flux equation q = —n c V T, 
where k c = nksKdynamic/electri n ls electron number density, and 
kg is the Boltzmann constant, for both electron and advective heat 
transport. 



(A/ijMJ 1 . For 1 < M A < (£/A) 1/3 the mean free 
path is less than I a which results in n e i ectr being some 
fraction of K, spitzer , while Kdynamic is given by Eq. ©. 
Thus Kdynamic/ Keiectr ~ j3aM s (L/\), i.e. the ratio does 
not depend on Ma (horisontal line in Fig. 1). When 
M A > (£A) 1/3 the mean free path of electrons is con- 
strained by l A . In this Case Kdynamic/ K e lectr ~ (3aM s M\ 
(see Eqs. and |@J) . This results in the separation 
line (3aM s ~ AfJ 3 in Fig. 1. 

4. TURBULENCE AND HEAT TRANSFER IN ICM 

It is generally believed that ICM is turbulent. The con- 
siderations below can be used as guidance. In unmagna- 
tized plasma with the ICM temperatures T ~ 10 s K and 
and density 10 -3 cm -3 the diffusivity vb=q ~ Vi on Xi on , 
where Vi on and Xi on are the velocity of an ion and its 
mean free path, respectively, would make the Reynolds 
number Re = LVl/v of the order of 30. This is barely 
enough for the onset of turbulence. For the sake of 
simplicity we assume that ion mean free path coin- 
sides with the proton mean free path and both scale as 
A w 3T 2 nZ\ kpc, where the temperature T3 = kT/3 keV 
and n_3 = n/10 -3 cm~ 3 . This provides A of the order 
of 0.8-1 kpc for the ICM (see NM01). 

It is accepted, however, that magnetic fields decrease 
the diffusivity. Somewhat naively assuming the maximal 
scattering rate of an ion, i.e. scattering every orbit (the 
so-called Bohm diffusion limit) one gets the viscosity per- 
pendicular to magnetic field v±_ ~ Vi on fLar,ioni which is 
much smaller than vb=o 1 provided that the ion Larmor 
radius rf, ari j 0n -C Xion - For the parameters of the ICM 
this allows essentially invicid motions 6 of magnetic lines 
parallel to each other, e.g. Alfven motions. 

In spite of the substantial progress in understading of 
the ICM (see Enfilin, Vogt & Pfrommer 2005, henceforth 
EVP05, Enfilin & Vogt 2006, henceforth EV06 and refer- 
ences therein), the basic parameters of ICM turbulence 
are known within the factor of 3 at best. For instance, 
the estimates of injection velocity Vl varies in the litera- 
ture from 300 km/s to 10 3 km/s, while the injection scale 
L varies from 20 kpc to 200 kpc, depending whether the 
injection of energy by galaxy mergers or galaxy wakes 
is considered. EVP05 considers an illustrative model in 
which the magnetic field with the 10 //G fills 10% of the 
volume, while 90% of the volume is filled with the field 
of B ~ 1 /xG. Using the latter number and assuming 
Vl = 10 3 km/s, L = 100 kpc, and the density of the hot 
ICM is 10" 3 cm" 3 , one gets V A « 70 km/s, i.e. M A > 1. 
Using the numbers above, one gets Ia ~ 30 pc for the 
90% of the volume of the hot ICM, which is much less 
than Xion- The diffusivity of ICM plasma gets v = Vi on lA 
which for the parameters above provides Re ~ 2 x 10 3 , 
which is enough for driving superAlfvenic turbulence at 
the outer scale L. However, as I a increases as oc B 3 , 
Re gets around 50 for the field of 4 fiG, which is at the 

6 A regular magnetic field B\ (2mkT) 1 ^ 2 c/(e\) that makes 
r Lar,ion l ess than A and therefore u± < vb=o is just 10 — 20 G. 
Turbulent magnetic field with many reversals over ri, ar i on does 
not interact efficiently with a proton, however. As the result, the 
protons are not constrained until I a gets of the order of r^ ar i on . 
This happens when the turbulent magnetic field is of the order of 
2 X 10 J (VL/10 3 km/s) G. At this point, the step for the random 
walk is ~ 2 X 10~ 6 pc and the Reynolds number is 5 X 10 10 . 
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border line of exciting turbulence 7 . However, the regions 
with higher magnetic fields (e.g. 10 //G) can support 
Alfvenic-type turbulence with the injection scale I a and 
the injection velocities resulting from large-scale shear 
V L (l A /L) ~ V L MX 3 . 

For the regions of B ~ 1 /iG the value of I a is smaller 
than the mean free path of electrons A. According to 
Eq. the value of Keiectr is 100 times smaller than 
Kspitzer- On the contrary, ^dynamic for the ICM param- 
eters adopted will be ~ "SQuspitzer, which makes the dy- 
namic diffusivity the dominant process. This agrees well 
with the observations in Voigt & Fabian (2004). Fig. 1 
shows the dominance of advective heat transfer for the 
parameters of the cool core of Hydra A ( B = 6 /iG, 
n = 0.056 cm -3 , L = 40 kpc, T = 2.7 keV according 
to EV06), point "CC", and for the illustrative model in 
EVP05, point "I", for which B = 1 /j,G. 

Note that our stationary model of MHD turbulence in 
§2 is not directly applicable to transient wakes behind 
galaxies. The ratio of the damping times of the hydro 
turbulence and the time of straightening of the magnetic 
field lines is ~ M^ 1 - Thus, for Ma > 1, the magnetic 
field at scales larger than I a will be straightening gradu- 
ally after the hydro turbulence has faded away over time 
L/Vl- The process can be characterized as injection of 
turbulence at velocity Va but at scales that increase lin- 
early with time, i.e. as lA+V A t. The study of heat trans- 
fer in transient turbulence and magnetic field "regularly" 
stretched by passing galaxies will be provided elsewhere. 

5. DISCUSSION AND SUMMARY 

In the paper above we attempted to describe the heat 
transfer by electron and turbulent motions for Ma < 1 
and Ma > 1. Unlike earlier papers, we find that turbu- 
lence may both enhance heat conduction and suppress it. 
For instance, when A gets larger than Ia the conductiv- 
ity of the medium ~ M^ 3 and therefore the turbulence 
inhibits heat transfer, provided that K e i ectr > Kdynamic- 
Along with the plasma effects that we mention below, 
this effect can, indeed, support sharp temperature gradi- 
ents in hot plasmas with weak magnetic field. 

As discussed above, rarefied plasma, e.g. ICM plasma, 
has large viscosity for motions parallel to magnetic field 

7 One can imagine dynamo action in which superAlfvenic tur- 
bulence generates magnetic field till Ia gets large enough to shut 



and marginal viscosity for motions that induce perpen- 
dicular mixing. Thus fast dissipation of sound waves 
in the ICM does not contradict the medium being tur- 
bulent. The later may be important for the heating of 
central regions of clusters caused by the AGN feedback 
(see Churasov et al. 2001, Nusser, Silk & Babul 2006 
and more references in EV06). Note, that models that 
include both heat transfer from the outer hot regions and 
an additional heating from the AGN feedback look rather 
promissing (see Ruszkowkski & Begelman 2002, Piffaretti 
& Kaastra 2006). We predict that the viscosity for 1 /iG 
regions is less than for 10 /xG regions and therefore heat- 
ing by sound waves (see Fabian et al. 2005) could be 
more efficient for the latter. Note, that the plasma in- 
stabilities in collisionless magnetized ICM arising from 
compressive motions (see Schekochihin & Cowley 2006, 
Lazarian & Beresnyak 2006) can resonantly scatter elec- 
trons and protons and decrease A for both species com- 
pared to the classical plasma values (A gets different for 
electrons and protons in this case). This decreases fur- 
ther K e i ec t r compared to n sp itzer but increases Re. In 
addition, we disregarded mirror effects that can reflect 
electrons back (see Malyshkin & Kulsrud 2001 and ref- 
erences therein), which can further decrease K e i ec tr- 

All in all, we have shown that it is impossible to charac- 
terize the heat transfer of magnetized plasma by a single 
fraction of Spitzer's value. The actual heat transport de- 
pends on sonic and Alfven Mach numbers of turbulence 
and may be much higher and much lower than the classi- 
cal one. As the result, turbulence can inhibit or enhance 
heat conductivity depending on the plasma magnetiza- 
tion and turbulence driving. Our study indicates that in 
many cases related to ICM the advective heat transport 
by dynamic turbulent eddies dominates thermal conduc- 
tivity. 
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